Trial wave functions with long-range Coulomb correlations for 
two-dimensional A^-electron systems in high magnetic fields 
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A new class of analytic wave functions is derived for two 
dimensional A'^-electron {2 < N < oo) systems in high mag- 
netic fields. These functions are constructed through breaking 
(at the Hartree-Fock level) and subsequent restoration (via 
post-Hartree-Fock methods) of the circular symmetry. They 
are suitable for describing long-range Coulomb correlations, 
while the Laughlin and composite-fermion functions describe 
Jastrow correlations associated with a short-range repulsion. 
Underlying our approach is a coUectively-rotating-electron- 
molecule picture, yielding for all A'' an oscillatory radial elec- 
tron density that extends throughout the system. 

Pacs Numbers: 73. 21. La, 73.43.-f, 73.22.Gk 



I. INTRODUCTION 

Two-dimensional (2D) few-electron systems in strong 
magnetic fields have been the focus of extensive theoret- 
ical investigations in the last twenty years. Many of 
these studies have used the Jastrow-Laughlin^ (JL) and 
composite-fermion'^ (CF) wave functions, where the dy- 
namics of electrons in extended fractional quantum Hall 
(FQH) systems is governed by the so-called Jastrow cor- 
relations, ft was shown^^ that the JL functions are ex- 
act eigenstates of the iV-electron problem under high- 
magnetic fields for a special short-range interparticle re- 
pulsion. However, based on close-to-unity overlaps with 
exact numerical solutions"'^'^"^'^ of the Coulomb problem 
for few-electron systems (with iV < 8), it is believed^^*'^^ 
that the JL/CF functions should not differ significantly 
from the exact Coulombic solutions. 

Recent experiments^^ on electron tunneling into the 
edges of a FQH system have found a current-voltage 
power law behavior, / cx V", with values for the expo- 
nent a that are in conflict with the universal prediction 
a = l/v derived from Jastrow correlations. These find- 
ings motivated^^'^* detailed exact diagonalization studies 
of FQH systems at the filling factor i' — 1/3 with up to 
A'' = 12 electrons. These latest studies revealed that the 
long-range Coulomb correlations lead to the formation 
of stripe-like oscillations in the radial electron densities 
(ED's) which are responsible^^ for the observed unex- 
pected behavior of the current-voltage power law. Most 
importantly, the JL functions fail^* to capture these ED 
oscillations, in spite of having overlaps with the exact 
wave functions that are very close to unity. 

For the iV-electron problem in strong magnetic fields 



and in the disk geometry (case of quantum dots, QD's), 
we use in this paper a microscopic many-body approach 
to derive analytic wave functions that capture the long- 
range correlations of the Coulomb repulsion. To obtain 
analytic results, we specifically consider the limit when 
the confining potential can be neglected compared to the 
confinement induced by the magnetic field. 

Underlying our approach is a physical picture of a col- 
lectively rotating electron molecule (REM) and the syn- 
thesis of the states of the system consists of two steps: 
First the breaking of the rotational symmetry at the 
level of the single-determinantal unrestricted Hartree- 
Fock (UHF) approximation yields states representing 
electron molecules (EM's, or finite crystallites). Sub- 
sequently the rotation of the electron molecule is de- 
scribed through restoration of the circular symmetry via 
post Hartree-Fock methods, and in particular Projec- 
tion Techniques'^ (PT's). Naturally, the restoration of 
symmetry goes beyond the mean-field and yields multi- 
determinantal wave functions. In contrast to the JL/CF 
functions, our analytic functions (applicable for any N 
and fractional filling) yield oscillatory radial ED's in 
agreement with the exact solutions of the A''-electron 
Coulombic system (see below). 



II. METHOD AND TRIAL WAVE FUNCTIONS 

In general, the symmetry-broken UHF orbitals are de- 
termined numerically. "'^'^'''''^°'^' However, in the case of 
an infinite 2D electron gas in strong magnetic fields, it 
has been found^^ that such UHF orbitals can be approx- 
imated by analytic Gaussian functions centered at differ- 
ent positions Zj = Xj + lYj and forming an hexagonal 
Wigner crystal (each Gaussian representing a localized 
electron). Such displaced Gaussians are written as (here 
and in the following i = V— T) 



l{z,Z,) = {1/^) 

X exp[-|z - Zjp/2] exp[-i(a;Y,- + yX^)] 



(1) 



where the phase factor is due to the gauge invariance. 
z = X — ly, and all lengths are in dimensionless units of 
IbV^ with the magnetic length being Is = y^hc/eB. 

In the case of a Coulombic finite iV-electron system, it 
has been found'^'^'^ that the UHF orbitals arrange them- 
selves in concentric rings forming EM's (referred to also 
as Wigner molecules, WM's).^'^ The UHF results for the 
ring arrangements are in agreement with the molecular 
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structures obtained via the conditional probability dis- 
tributions (CPD's, which can be extracted from exact 
numerical wave functions^"'^^'^^), as well as with those 
obtained^^ for the equilibrium configurations of classical 
point charges in a 2D harmonic trap.^^ 

For an TV-particle system, the electrons are situated at 
the apexes of r concentric regular polygons. The ensu- 
ing multi-ring structure is denoted by (ni, n2, n^) with 
X]o=i ~ ^- The position of the j-th electron on the 
q-th ring is given by 



Z« = Zg exp[i27r(l - i)/n,], 1 < j < rig 



(2) 



We expand now the displaced Gaussian (1) over the 
Darwin- Fock single-particle states. Due to the high mag- 
netic field, only the single-particle states. 



exp(— 22;*/2) , 



(3) 



of the lowest Landau level are needed (the angular mo- 
mentum of this state is — / due to the definition z = x—iy) 
Then a straightforward calculation^^ yields 



u{z,Z) = Y,Ci{Z)Uz) , 

1=0 



(4) 



with Ci{Z) = {Z*yexp{-ZZ*/2)/VV. for Z ^ 0. Natu- 
rally, Co(0) = 1 and Ci>o(0) = 0. 

Since electrons in strong magnetic fields are fully po- 
larized, only the space part of the many-body wave func- 
tions needs to be considered.^* The symmetry-broken 



UHF determinant, 4" 



N 

UHF) 



describing the WM is con- 



structed out of the localized wave functions u{z, Zj). Us- 
ing (4) one finds the following expansion (within a pro- 
portionality constant) 



*UHF 



E 

Ii=0,...,In=0 



Ci^{Zi)Ci^{Z2) ■ ■ ■ Ci^{Zm) 

N 



X D{h, h, In) exp(- ^ ZiZ^/2) 



(5) 



i=l 



where D{li,l2, --Jn) = dct[z[\ z!^^. ■ ■ ■,z'j^]. The Zk's 
(with l<k<N)m Eq. (5) are the Z]'s of Eq. (2), but 
relabeled. 

The UHF determinant ^'uhf breaks the rotational 
symmetry and thus it is is not an eigenstate of the total 
angular momentum JiL = h'^f^-^li. However, one can 
restore^^'^^ the rotational symmetry by applying onto 
^UHF following projection operator^^ 



/■27r 



where hL. 

,9-1 



djg eXp[l'yg{Lg - Lg)] , (7) 

^^lii and hLq = with 

lg = ^l-i Tig {ii — 0) are partial angular momenta op- 
erators and values, respectively, associated with the gth 
ring, and TiL = h J2g^i Lg are the eigenvalues of the total 
angular momentum. 

When applied onto ^'uhf' ^^'^ projection operator Ol 
acts as a product of Kronecker deltas: from the unre- 
stricted sum (5), it picks up only those terms having a 
given total angular momentum L and a specific ordered 
partition of it into partial angular momenta associated 
with the concentric rings, i.e., UL = ti Y^^g^i Lq. The final 
analytic expression depends on the specific ring arrange- 
ment (ni, 712, rir). For lack of space, we will present 
here explicitly only the simplest nontrivial arrangement, 
i.e., (711,712), with more complex [or simpler ones, i.e., 
(0,N) and (1, N ~ 1)] obtained via straightforward exten- 
sions. 

For specific electron locations (2) associated with the 
(711,712) WM, one derives^'' the following symmetry- 
preserving, many-body correlated wave functions (within 
a proportionality constant), 

*z,i,z,2('^i>"-2; [z]) = 

iiH l-i„i=ii, i„,i + iH \-In=L2 /n 

0<h<l2<---<lN 



XI JJ^ sin 

. l<i<j<ni 



XI sin 

yni + l<i<j<N 



?^2 



N 



X 2?(Zi,«2,...,Mexp(-X]^i4/2) • (8) 

i=l 

In deriving (8), we took into account that for each de- 
terminant D{li,l2, ■■■,In) in the unrestricted expansion 
(5) there are A'^! — 1 other determinants generated from 
it through a permutation of the indices {Zi, ^2) ^at}; 
these determinants arc equal to the original one or differ 
from it by a sign only. 

Generalizations of expression (8) to structures with a 
larger number r of rings involve for each additional q-th 
ring {2 < q < r): (I) Inclusion of an additional product 
of sines with arguments containing 77^; (II) A restriction 
on the summation of the associated angular momenta. 
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with 



III. PROPERTIES OF THE REM WAVE 
(6) FUNCTIONS 

We call the correlated wave functions [Eq. (8)] the 
REM wave functions. Among the properties of the REM 



2 



TABLE I. The Ql[z] polynomial associated with the 
EMWF's and the JL functions (The [z] polynomials are 
of order L — Lo). 



EMWF 




- 3ZiZ2 + z'i + 6Z1Z2Z3 


- 32I23 - 


32i2| + zi) 








- 3ziZ2 + z'i + 6212223 


— 32i23 - 


3z2zl + zi) 


1 DD 


JL 




{zi - 22)^(21 - 23 


"(22-23^ 







functions, we mention the following: 

1) The REM wave functions lie entirely within the 
Hilbert subspace spanned by the lowest Landau level and, 
after expanding the determinants,'^'^ they can be written 
in the form (within a proportionality constant), 



N 



$^[z] = Pf[z]exp(-5]z,z,r/2), 



(9) 



i=l 



where the [z] 's are order-L homogeneous polynomials 
of the Zi's. 

2) The polynomials Pl'iz] are divisible by 



n (-^ 

l<i<j<N 



(10) 



namely Pj^ [z] = Py [z]Qy^[z]. This is a consequence of 
the antisymmetry of <I>;^ [z] . Py [z\ is the Vandermonde 
determinant I?(0, 1, N). For the case of the lowest al- 
lowed angular momentum Lq = N{N — l)/2 (see below), 
one has P/^ [z] = Py [z] , a property that is shared with 
the Jastrow-Laughlin^ and composite-fermion^ trial wave 
functions. 

3) The P^[zys are translationally invariant functions. 

4) The coefficients of the determinants [i.e., products 
of sine functions, see Eq. (8)] dictate that the REM func- 
tions are nonzero only for special values of the total angu- 
lar momentum L given for a (ni, n2, n^) configuration 
by, 



N{N - l)/2- 
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0,1,2,3,... (11) 



The minimum angular momentum Lq = N(N — l)/2 is 
determined by the fact that the D determinants [see Eq. 
(8)] vanish if any two of the single-particle angular mo- 
menta k and Ij are equal. For the (0, N) and (1, iV — 1) 
rings, the special values are given hy L — Lq + Nk and 
L = Lo + {N—l)k, respectively. In plots of the energy vs. 
the angular momenta, derived from exact-diagonalization 
studies, ^""'^^ it has been found that the special L values as- 
sociated with the (0, N) and (1, N— 1) rings (appropriate 
for N < 7) exhibit prominent cusps reflecting enhanced 
stability; as a result these L values are often referred to 
as "magic angular momenta" .^^ We predict that similar 
magic behavior reflecting enhanced stability is exhibited 
by the special L values given by Eq. (11) and associ- 
ated with the general ring arrangement (rii, n2, n^). 
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FIG. 1. Radial ED's from EMWF's (solid lines, all frames), 
exact diagonalization [dashed lines, (a) — (e)], and JL functions 
[dotted lines, (d) and (e)]. In (a) and (b), the solid and dashed 
curves are practically indistinguishable. 



In the thermodynamic limit, the total L is related to 
a fractional filling f — N{N — l)/(2i), and thus the an- 
gular momenta (11) of the REM functions correspond to 
all the associated with the FQHE, including the even- 
denominator ones, i.e., v = I, 3/5, 3/7, 5/7, 2/3, 1/2, 
1/3, etc... 

5) For the case of two electrons {N — 2), the REM 
functions reduce to the Jastrow-Laughlin form, namely 



n 

i<i<j<Ar 



(12) 



where L — 1, 3, 5, ... However, this is the only case 
for which there is coincidence between the REM and 
the JL wave functions. For higher numbers of electrons, 
N , the polynomials P^[z] of the REM functions (apart 



from the lowest-order Vandermonde [z] ones) are 
quite different from the corresponding JL or composite- 
fermion polynomials. In particular, the familiar factor 
Y\i<^i^j<,jq{zi — Zj)'^P , withp an integer,^-'* (which reflects 
multiple zeroes) does not appear in the REM functions 
(see, e.g.. Table I which contrasts the Qg[z] polynomials 
corresponding to the REM and JL functions). 

6) For the case of three electrons {N = 3), after trans- 
forming to the Jacobi coordinates z = (zi + Z2 + z^)/^, 

Za = (2/3)l/2((zi + Z2)/2 - Z3), Zi, = (Zi - Z2)lV2 

(and dropping the center-of-mass exponential factor), the 
REM wave functions can be written as (again within a 
proportionality constant), 

$|[Za, Zh] = [{Za + IZh)^ ~ {Za ~ IZb)^] 

xexp[(-l/2)(za< + ^fc4)] , (13) 

with L = 3m, m = 1, 2, 3, 4, ... being the total an- 
gular momentum. Again the wave functions $|^[zQ,Zb] 
are very different from the three-electron JL ones; e.g., 
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they arc nonvanishing for even m values, unlike the three- 
electron JL functions. However, the ^^[za, ztYs coincide 
with the functions |m, 0) derived in Ref. 1. We remark 
that, although it was found^ that these wave functions 
exhibited behavior expected of fractional quantum Hall 
ground states (e.g., areal quantization and incompress- 
ibility), the generalization of them to a higher number 
of electrons did not follow; instead, the REM functions 
presented here do constitute such a generalization. 

IV. OSCILLATORY ELECTRON DENSITIES 

In Fig. 1, we display the radial ED's of several REM 

wave functions and compare them to corresponding ED's 
from exact diagonalizations. The main conclusion is that 
the ED's of the REM functions exhibit^^ a prominent os- 
cillatory behavior in excellent agreement with the exact 
ED's. Such an oscillatory behavior is a natural conse- 
quence of the underlying ring arrangements. For N = 6 
and L = 20, the underlying structure is a (1,5) arrange- 
ment [L is 5 units larger than the minimum Lq = 15, i.e., 
ni = 1, fci = and n2 = 5, fc2 = 1 in Eq. (11)], and thus 
the corresponding ED exhibits a maximum at the origin 
followed by an outer hump [Fig. 1(a)]. For N = 6 and 
L = 21, however, the underlying structure is a (0,6) ar- 
rangement {L is 6 units larger than Lq), and thus the ED 
exhibits a dip at the origin and a single outer hump [Fig. 
1(b)]. In the other = 6 cases plotted here [Figs. 1(c), 
1(d), 1(e)], the difference L — Lq is divisible by 5 and the 
underlying ring arrangement is (1,5); thus there are two 
humps in the corresponding ED's, the inner one portray- 
ing the single electron at the origin. The N = 6, L = 75 
[Fig. 1(d) ] and iV = 6, L = 105 [Fig. 1(c)] cases corre- 
spond to the fractional fillings 1/5 and 1/7, respectively. 
For these two cases, we have plotted also the ED's associ- 
ated with the JL fimctions (dotted lines). As was found 
in the u = 1/3 case,^^ the JL functions fail to capture 
the radial oscillations that are characteristic of the long- 
range Coulomb force. Finally, in Fig. 1(f) and Fig. 1(g), 
we present the ED's of the REM functions for N = 12, 
L = 123 and N = 16, L = 205. In general, there are 
as many humps as the number of concentric rings. In- 
deed the ring structures are^^'^^ (3,9) and (1,5,10) for 
= 12 and A^ = 16, respectively. ■'^■^ We note that the 
ED for N = 12, L = 123 is similar to the exact ED for 
the same case,^^ although the latter was calculated with 
an external confinement. 

It has been found^^''^^ that the number of rings (r) in- 
creases as the number of electrons grows. For v < 1/3, 
this results in electron-density oscillations that extend 
along the whole radius of the QD, with no obvious sepa- 
ration into bulk and edge regions. Currently, the largest 
number of electrons for which the ring structure has been 
determined^^ is N = 230 with a concentric-ring arrange- 
ment of (1,6,12,18,23,25,34,37,37,37). 



V. CONCLUSIONS 

We have developed a new class of trial wave functions 
of simple functional form, which accurately describe the 
physics of electrons in QD's under high magnetic fields. 
In particular, our functions capture the long-range cor- 
relations of the Coulomb repulsion; unlike the JL func- 
tions, they yield for all A^ and fractional fillings an os- 
cillatory radial electron density in agreement with exact- 
diagonalization results. The electron density oscillations 
extend throughout the system. The thematic basis of 
our approach is built upon the intuitive, but microscop- 
ically supported, picture of collectively rotating electron 
molecules, and the synthesis of the many-body REM 
wave functions involves breaking of the circular symme- 
try at the UHF level with subsequent restoration of this 
symmetry via a projection technique. While we focus 
here on the strong magnetic-field regime, we note that 
the REM picture unifies the treatment of strongly corre- 
lated states of electrons in QD's over the whole magnetic- 
field range. Finally, our REM wave functions, 
aimed here mainly at treating finite electron systems (i.e., 
QD's), can provide in the thermodynamic limit an al- 
ternative interpretation of the FQH effect; namely, the 
observed hierarchy of fractional filling factors may be 
viewed as a signature originating from the magic angular 
momenta of rotating electron molecules. 

This research is supported by the U.S. D.O.E. (Grant 
No. FG05-86ER-45234). 
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